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ABSTRACT 

I summarise recent lattice results on the QCD confinement mechanism in the 
maximally Abelian projection. 



1. Introduction 

The phenomenology of strong interactions contains two fundamental ingredients: 
asymptotic freedom and the confinement of colour charges. The former requirement led 
to the invention of QCD. A mathematically rigorous proof that QCD as the microscopic 
theory of strong interactions indeed gives rise to the macroscopic property of linear 
quark confinement as indicated by Regge trajectories and quarkonia spectra is, after 
a quarter of a century, still lacking. Meanwhile, lattice gauge theory simulations have 
provided convincing numerical evidence for this conjecture. 

The difficulty in deriving infra red properties of QCD illustrates that something 
qualitatively new is happening: unlike in previously existing elementary physical theo- 
ries, it is not possible to reduce everything down to two-body interactions but collective 
excitations of quark and gluon states have to be accounted for. For the first time, ex- 
citations of the vacuum that are considered to be fundamental do not occur as initial 
or final states anymore. Therefore, understanding confinement, in my opinion, is one of 
the most exciting challenges of modern physics. New physical and mathematical tech- 
niques that can successfully be applied to non-perturbative QCD might be required 
for dealing with other strongly interacting theories or theories with a non-trivial vac- 
uum structure, in general. Vice versa techniques developed in a different context might 
help to prove QCD confinement and in solving QCD. A recent example is the proof 
of confinement in SUSY Yang-Mills theories.0 Of course QCD as we know it does not 
obey super-symmetry but nonetheless such activities point into a promising direction. 
Just as QCD can serve as a guinea pig and development centre for non-perturbative 
techniques, lattice simulations help probing the validity range of effective low-energy 
models or in verifying certain conjectures. 

QCD in itself is sufficiently complicated to keep many physicists busy. Solving 
QCD is still important and even more so since (almost) everyone believes in it. It 
is widely accepted that perturbative QCD (pQCD) successfully describes high en- 
ergy scattering processes. However, without understanding non-perturbative aspects 
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of QCD, even in this case, it is not clear why pQCD works at all. Only hadrons, lep- 
tons and photons but no quarks or gluons are observed in experiment. In order to 
fill the gap, the formation of colour-neutral hadronic jets from quarks and gluons has 
to be modelled. Furthermore, it is commonly conjectured that, once the fundamental 
scattering on the quark and gluon level has taken place, further interactions can be 
neglected. It is demanding to derive the ingredients of fragmentation models directly 
from truly non-perturbative QCD, as well as to verify the factorisation hypothesis. 

Many phenomenologically important questions are posed in low-energy QCD that 
eagerly await an answer: is the same set of fundamental parameters (QCD coupling and 
quark masses) that describes for instance the hadron spectrum consistent with high 
energy QCD or is there place for new physics? Are all hadronic states correctly classified 
by the naive quark model or do glueballs, hybrid states and molecules play a role? At 
what temperatures/densities does the transition to a quark-gluon plasma occur? What 
are the experimental signatures of quark-gluon matter? Can we solve nuclear physics 
on the quark and gluon level? Clearly, complex systems like iron nuclei are unlikely 
ever to be solved from first principles alone but modelling and certain approximations 
will always be required. Therefore, it is desirable to test model assumptions, to gain 
control over approximations and, eventually, to derive low-energy effective Lagrangians 
from QCD. Lattice simulations are a very promising tool for this purpose. 

In the past decades, many explanations of the confinement mechanism have been 
proposed, most of which share the feature that topological excitations of the vacuum 
play a majoiLrple. A list of these theories, includes the dual superconductor picture of 
confinement ,0a the centre vortex model,a the instanton liquid model,0 and the anti- 
ferromagnetic vacuum.Q All these interpretations have been explored in lattice studies. 
The situation with respect to an anti-ferromagnetic vacuum is still somewhat incon- 
clusiveJ3 Instantons seem to be more vital for chiral symmetry related properties than 
for confinement Depending on the picture, the excitations giving rise to confinement 
are thought to be magnetic monopoles, instantons, dyons, centre vortices, etc.. I would 
like to stress that the above ideas are not completely disjoint and do not necessarily 
exclude each other. For instance, all the above mentioned topological excitations are 
found to be correlated with each other in numerical as well as analytical studies. 

I will restrict myself to the, at present, most popular superconductor picture which 
is based on the concept of electro-magnetic duality after an Abelian gauge projection. 
Recently, the centre vortex model has been resuscitated too. In the latter picture, 
the centre group that is directly related to the traditional order parameter of the de- 
confinement phase transition in finite temperature pure gauge theories, the Polyakov 
line, plays an essential role. Another striking feature is that — unlike monopole cur- 
rents — centre vortices form two-dimensional objects, such that in four space-time 
dimensions, a linking number between a Wegner- Wilson loop and centre vortices can 
unambiguously be defined, providing a geometric interpretation of the confinement 
mechanism.a Unfortunately, I have not enough space to review these exciting develop- 
ments. Therefore, I refer the interested reader to a review by Greensite.cJ 

Another obvious omission are questions related to chiral properties of QCD. Real 
world QCD is not only a theory of gluons but includes quarks which means that the de- 
confinement phase transition will eventually be replaced by a chiral phase transition. In 



the phase with broken chiral symmetry colour charges are still being anti-screened and 
linear confinement approximately holds. However, if the binding energy within a hadron 
exceeds a critical value, the hadron will break up into two or more colour-neutral parts 
(string breaking). 

Based on different confinement pictures various effective models of infra red QLTD 
have been proposed in the past. Three such examples are the Abelian Higgs modeljiim 
the stochastic vacuum model£3 and dual QCD HI As well as understanding the mecha- 
nism of confinement, it is desirable to verify aspects of these models and eventually to 
derive certain model parameters directly from QCD. 

This article is organised as follows. I will start with a brief motivation and in- 
troduction into lattice methods in Sec. |2|, before explaining the concept of the Abelian 
projection in Sec. ||. Subsequently, I will introduce some notations and review lattice 
results in Sec. |] and conclude with a summary of answered and open questions. 

2. Why Lattice? 

Lattice methods allow for a somewhat brute force but first principles numerical 
evaluation of expectation values of a given quantum field theory that is defined by an 
action S. The lattice volume and spacing are limited due to finite computer speed and 
memory. Simulations are performed in Euclidean space and an analytic continuation 
to Minkowski space of numerical results that have been obtained on discrete points 
with finite precision is virtually impossible. Therefore, not every physically meaningful 
quantity can be calculated in a straight forward manner. The obvious strength of lattice 
methods are hadron mass determinations. Even if one is unimpressed by post-dictions 
of experimental values that have been known with high precision for decades, within 
an accuracy of only 5 %, such simulations allow one to fix fundamental standard model 
parameters like quark masses and QCD running coupling in the low-energy domain. 
Of course, plenty of other applications of phenomenological importance exist. 

Unfortunately, only the lowest lying one or two radial excitations of a hadronic 
state are accessible in practice. Lattice predictions are restricted to rather simple sys- 
tems too. Even the deuterium is beyond the reach of present day super-computers. 
Therefore, it is desirable to supplement lattice simulations by analytical methods. The 
computer alone acts as a black box. In order to understand and interprete the output 
values and to predict their dependence on the input parameters, some modelling is 
required. Vice versa, the lattice itself is a strong tool to validate models and approx- 
imations. Unlike in the "real" world, we can vary the quark masses m q , the number 
of colours N c , the number of flavours rif, the temperature, the volume, the space-time 
dimension of our /emto-universe and even the boundary conditions in order to expose 
models to thorough tests in many situations. Instead of indirectly and in a somewhat 
uncontrolled fashion deriving parameter values from experiment we can compute cus- 
tom designed observables. What is thought to work in the real world ought to work on 
the lattice too! Moreover, many models rely on certain approximations. Experimental- 
ists cannot switch off quark flavours but we can! 

In a lattice simulation, hyper-cubic Euclidean space-time is discretised on a box 
with, say, L 4 lattice points or sites, x. Two adjacent points are connected by an oriented 
bond or link (x,fi). The product of four links, enclosing an elementary square, is a 



plaquette. Quarks are living on sites, gauge fields on links and the plaquette determines 
the curvature within the SU(N C ) group manifold and corresponds to the field strength 
tensor. For simplicity, we assume an isotropic lattice with equal lattice spacing a in all 
directions. The lattice spacing provides an ultra violet cut-off on the gluon momenta 
q < it /a and regulates the theory. 

We simulate an action S(/3, m l q ) on the lattice, which contains the quark masses 
as well as an inverse QCD coupling (3 = 2N C / g 2 as free parameters. By varying j3 (and 
m l ) the lattice spacing a is changed. Asymptotic freedom tells us that the ultra violet 
cut-off is removed as (3 — > oo. The physical dimension of a is determined by calculating 
a dimensionful quantity on the lattice and associating it to its experimental value. If the 
right theory is being simulated we should be able to reproduce all experimental mass 
ratios in the continuum limit a — > 0, such that it becomes irrelevant what experimental 
quantity we have chosen to set the scaleQ. In practice one does not get all the quark 
masses right, such that there are always systematic uncertainties in a, due to the 
ambiguity of the choice of this experimental input quantity. 

The discretised lattice action is formulated in a manifestly gauge-invariant way 
and approaches the continuum action with a —>■ 0. One ideally extrapolates to this limit 
at fixed physical lattice extent, La = const.. Later on, the thermodynamic limit La —>■ 
oo should be investigated. Due to the finiteness of a computer these two extrapolations, 
as well as extrapolations and interpolations between results obtained at different quark 
masses, are subject to systematic uncertainties that have to be carefully estimated. 
For most hadronic processes, a lattice spacing a < 0.1 fm is considered to be close to 
the continuum limit while an extent La ~ 2 fm is comfortably large to accommodate 
ground state mesons and baryons. 

Expectation values of operators O are determined by computation of the path 
integral, 

{0) = \ I [ df/ H#][#] 0[U]e- s ^. (1) 

The normalisation factor Z is such that (1) = 1. U Xjtt G SU(N C ) denotes a gauge 
field and the 4 x n/ x N c tuple ip x is a fermion field. The high- dimensional integral 
is evaluated by means of a (stochastic) Monte-Carlo method as an average over an 
ensemble of n representative gauge configurations, Q = {U®}, i — 1, . . . , n. Therefore, 
the result on the expectation value is subject to a statistical error that will basically 
decrease like 1 / y/n: the longer we measure the more precise the prediction becomes. 
Therefore, we might speak of lattice measurements and lattice experiments in analogy 
to "real" experiments. Ideally, the sample size n is such that the statistical precision is 
similar or smaller than the systematic uncertainty of the extrapolations. 

In Fig. [I] I illustrate that confinement is a numerical fact by displaying the poten- 
tial between a static quark and an anti-quark, separated by a distance r. The data have 
been obtained in 577(3) gauge theory with and without two light (mass-degenerate) 

*In the results presented here, I set the scale (somewhat arbitrarily) by the value ^/a = 440 MeV for 
the string tension in case of SU(2) and by the value ro = 0.5 fm in case of SU(3) chromodynamics. ro 
is the distance at which dV{r)/dr\ rQ — 1.65 where V(r) denotes the potential between two static 
colour sources, separated by a distance r. 
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Fig. 1. Ground state (full symbols land E u hybrid potential (open symbols) at /3 = 6.2 
(quenched) and (3 = 5.6, K = 0.1575.E-S 

quark flavours at lattice spacings somewhat smaller than 0.1 fm. The scale corre- 
sponds to 0.5 fm. The (partially) unquenched potential, that has been obtained at a 
quark mass m q ~ 50 MeV, exhibits a more pronounced singularity at small r, due to the 
slower running of the QCD coupling towards zero as the momentum scale is increased 
(or r is decreased). Furthermore, from r « 1.2 fm onwards, the (partially) unquenched 
potential eventually flattens: the QCD string "breaks". Indeed, many interesting fea- 
tures of the theory can be studied on the lattice that are not at all accessible to real 
experiment^As an example, I have included a so-called E u (or continuum U u ) hybrid 
potentialtSlij into the plot that corresponds to the interaction energy of two static 
quark sources for the case that gluons contribute one unit of angular momentum along 
the interquark axis. 

3. (Abelian) photons and monopoles 

Soon after the advent of QCD, 't Hooft and MandelstanJ! proposed the dual 
superconductor scenario of confinement; the QCD vacuum is thought to behave analo- 
gously to an electrodynamic superconductor but with the roles of electric and magnetic 
fields being interchanged: a condensate of magnetic monopoles expels electric fields from 
the vacuum. If one now puts electric charge and anti-charge into this medium, the elec- 
tric flux that forms between them will be squeezed into a thin, eventually string-like, 
Abrikosov-Nielsen-Oleson (ANO) vortex which results in linear confinement. 

In all quantum field theories in which confinement has been proven, namely in 
compact U(l) gauge theory in the Villain formulationJIZi the Georgi-Glashow modeled 
and SUSY Yang-Mills theories,!^ this scenario is indeed realised. Before one can ap- 
ply this picture to QCD or SU(N C ) chromodynamics one has to identify the relevant 
dynamical variables: it is not straight forward to generalise the electro-magnetic du- 
ality of a £7(1) gauge theory to SU(N C ) since gluons carry colour charges. How can 



one define electric fields in a gauge invariant way? What fields are dual to the electric 
fields? How can one identify a monopole current? The next question is that for the 
order parameter of the confinement- deconfinement phase transition: what symmetry 
is broken? One might also ask why this symmetry is broken — a question that has 
been answered by the BCS theory for conventional superconductors. Finally, we are 
interested in an effective low-energy theory which corresponds to the Ginzburg-Landau 
theory of standard superconductors. 

Let us address the first question. In the Georgi-Glashow model, the SO (3) gauge 
symmetry is broken down to a U(l) symmetry as the vacuum expectation value of 
the Higgs field becomes finite. Within this effective U(l) theory, the standard electro- 
magnetic duality is realised, resulting in confinement. Such a mechanism is not provided 
by QCD but one can attempt to reduce the SU(N C ) symmetry to an Abelian symmetry 
by hand. In this spirit, it has been proposedQ to identify the monopoles in the {/(l)^ -1 
diagonal Cartan subgroup of SU(N C ) gauge theory after gauge fixing in respect to the 
off-diagonal SU(N c )/U(l) Nc ~ 1 degrees of freedom (Abelian projection). 

In general, a gauge transformation, Q x G SU(N c )/U(l) Nc ~ 1 , can be found that 
diagonalises an arbitrary operator X x = f[U ] within the adjoint representation of 
SU(N C ) when applied to the link variables, 

u% = ii r c,.,,n:..„ ir (2) 

Subsequently, a coset decomposition of the gauge transformed links is performed, 

^ = Ca with C x ^eSU(N c )/U(l) N °-\ Uv^eUil)"*- 1 . (3) 
If we now apply a residual gauge transformation, uo x G U{l) Nc ~ l , we find, 

Ux,fi — ► v x u x ,nul +a[l , C X)fl — ► L) x u x ,pul, (4) 

i.e. u x ^ transforms like a gauge field while C XifJi transform like matter fields. Therefore, 
we will refer to the diagonal gluon field as "photon" field. While the gauge fields 
of electrodynamics in the non-compact continuum formulation are free of singularities 
and, therefore, cannot form magnetic monopole solutions, the gauge transformation fl x 
is singular wherever two eigenvalues of X x coincide. Therefore, the gauge field u x ^ will 
in general contain such monopoles. 

After Abelian gauge fixing QCD can be regarded as a theory of interacting pho- 
tons, monopoles and matter fields (i.e. off-diagonal gluons and quarks). One might as- 
sume that the off-diagonal gluons do not affect long range interactions. This conjecture 
is known as Abelian dominances^ In addition, monopole dominance of non-perturbative 
physics has been proposed.^ 

The identification of photon fields and monopoles is a gauge invariant process. 
However, the choice of the operator X x , that defines the Abelian projection, is am- 
biguous. In his original work, 't Hooft suggested that the dual superconductor sce- 
nario would be realised in any Abelian projection. Indeed, the expectation value of 
a monopole creation operator has been found to be an order parameter of the de- 
confinement phase transition in quite a few different projections of SU(2) chromody- 
namics." As long as any ^(l)^ -1 projection of the theory yields similar results, one 



does not have to specify the mechanism that is thought to break the SU(N C ) gauge 
symmetry of QCD. However, numerical simulations suggest that Abelian and monopole 
dominance are not at all universal. The most popular ^auge projection applied and dis- 
cussed is the maximally Abelian projection (MAP). til One feature that is not shared 
by almost all other projections that have been investigated on the lattice so far is that 
the (local) MAP gauge condition gives rise to non-propagating ghost fields only, guar- 
anteeing renormalisability of the Abelian projected theory. Another — and possibly 
related — fact is that both, Abelian and monopole dominance have qualitatively been 
verified in this projection. 

One should mention that the analogy to an ordinary superconductor after Abelian 
projection is not complete. The electrons that form Cooper pairs in BCS theory are all 
negatively charged. However, the monopoles that are thought to condense in QCD can 
carry both, negative and positive magnetic charges. Therefore, the composition of the 
condensate is very different. The origin of the interaction that results in an attractive 
force between monopoles and anti-monopoles also differs from the periodic background 
potentials of the BCS theory. 

I will briefly explain how MAP is performed on the lattice for the case of SU (2) 
gauge theory. In the first step one maximises the functional,E3 

F(n)=J2tr(T 3 U^T 3 U?l), (5) 

by means of a gauge transformation VL X . T3 denotes a Pauli-matrixfj]. After the max- 
imisation, all link variables are as diagonal as possible. The resulting gauge fields 
Ux,^ = ex P {jl2cA c X41 T c /2 s j satisfy 't Hooft's differential MA gauge fixing condition,! 




After gauge fixing, the projection is performed: observables are calculated on 
Abelian configurations {0 X ^} rather than {U x ^}. We refer to the field theory, defined 
in this way, as Abelian projected SU{2) gauge theory [APSU(2)]. For convenience the 
Abelian links are represented in the Lie algebra rather than the group itself, u x ^ = 
exp {i9 x ^Ts). Note that the normalisation has been chosen such that the periodicity of 
Op, is 2tc as opposed to the periodicity 4-7T of A 3 ^. With this convention we find magnetic 
charges to be multiples of 2ir/e just like in electrodynamics rather than 2N c 7r/e as one 
might have expected in SU(N C ) gauge theory. We also use the convention e = 1 for the 
fundamental electric charge. 

Recently, some articles have appeared whose authors claim to have "proven" 
Abelian dominance, either in general or specifically in the MAP. In one case the argu- 
ment is based on the fact that the off-diagonal gluon fields acquire mass and, therefore, 

tin principle, we can fix the gauge along any direction within the SU{2) group space. The resulting 
gauges will only differ from each other by a global gauge transformation that will not affect expectation 
values — as long as we perform the subsequent coset decomposition with respect to the same U(l) 
subgroup. Also note that if we forced an adjoint Higgs field with a <5-like potential into the 3-direction, 
the form of the interaction term with the gauge fields would be identical to F . 



are thought to affect ultra violet physics only, a statement that does not necessarily ap- 
ply to confining field theories. Other arguments are either based on misinterpretations 
of the transfer matrix formalism or on rewriting the original theory in terms of Abelian 
projected variables plus perturbations (the off-diagonal fields) that are proportional 
to a small parameter e. In order to avoid going too much into details I only mention 
that both, the static potential and Wilson loops in MAP are very different from their 
counter parts in the original theory, even for large distances. In case of the potential, 
the asymptotic slope (string tension) comes out to be similar. However, APSU(2) and 
SU (2) potentials are shifted with respect to each other by a substantial constant due 
to different self-energy contributions. Any "proof" that fails to account for the latter 
fact is necessarily wrong to some extent. 

4. Results 

4-1. Some Definitions 

The language of differential formsil turns out to be very convenient for the present 
purpose and this is even more so on the lattice. Let us start with a U(l) field theory 
in D = 4 dimensions. The anti-symmetric field strength tensor F (a 2-form) can be 
decomposed as follows, 

F = dA + S*C, (7) 

where A denotes the standard magnetic four-potential (1-form) and C denotes an 
electric four-potential which is a 1-form on the dual lattice. A and C are not uniquely 
determined but subject to a £4/(1) X U mag {l) gauge invariance since gradients of scalar 
fields can be added, "d" denotes the exterior derivative that, applied to an n-form, 
results in an n + 1-form. "*" is the pull-back operator that connects an n-form to 
a (D — n)-form on the dual lattice while the dual derivative 5 = *d* turns n-forms 
into (n — l)-forms. We have * 2 = 1 and d 2 = 5 2 = 0. We define a Laplacian □ = 
(-) D (d5 + 5d). 

In this notation, we obtain the generalised Maxwell equations, 

5F = j, dF = *k. (8) 

The electric current j is a 1-form on the original lattice while the magnetic (monopole) 
current k is a 1-form on the dual lattice. In Landau gauge (5A = or 5C = 0, 
respectively) this means, 

j = DA, k = DC. (9) 

Note that, unlike the continuum four-potential A, the link angles 9 £ (— tt, tt] are 
subject to a 2tc shift-periodicity. Therefore, the identification, 

F= -4sin(d0) \l + 0(a 2 )] , (10) 
a 2 L J 

is natural for electro-magnetic fields. We can factorise the plaquette dO G (— 47r, 4ti} 
into a regular part 6 n £ (— tt, tt] and a singular part m G {—2, —1, 0, 1, 2}: 



d9 = 9 n + 2irm. 



(11) 



While ddd = 0, d9 n = —2iTdm does not necessarily vanish. Let us consider the magnetic 
flux through a (spatial) plaquette, $ ma5 = a 2 d9 [1 + 0(a 2 )}, which, in the absence 
of monopole strings, would be identical to a contour integral of the vector potential 
A around the plaquette. However, in presence of monopoles, only exp (i J n d 2 f d9) = 
exp (i J dn dx A) holds: the argument can be shifted by multiples of 2ir which, with 
the normalisation e = 1, is the elementary monopole charge. Thus, —2nm counts 
the monopole contribution to the magnetic flux through the plaquette and —2ndm 
corresponds to the "flux" out of a 3-cube, such thatEJ 

k = —2tc *dm = *d6 n (12) 

is a magnetic monopole current^, k is conserved since 5k = —2ix8 2 *m = 0, i.e. monopole 
currents form closed loops on the dual lattice. The individual components of k can take 
values 27m with n = —4, —3, . . . , 4. 

Monopole currents can alternatively be defined through,!! k = *dF = *dsin (d9), 
the advantage being that the second Maxwell equation [Eq. @] is automatically ful- 
filled. The current k is obviously conserved too. The quantisation of magnetic charges, 
however, is obscured by lattice artefacts. Magnetic monopoles become extended objects 
and a geometric interpretation is not as straight forward as for the definition presented 
before. For all results I am going to review, the first definition has been used. Note that 
the locations of monopoles in either of the definitions are not necessarily identical to 
positions of singularities of the adjoint operator X that has been diagonalised by the 
gauge fixing, as originally proposed by 't Hpoft. 

Each link variable can be factorisedcBca into a singular part stn9 that is induced 
by magnetic monopoles and a regular (or photon) part 9 reg = 9 — 9 stng that obeys 
the homogenous Maxwell equation dF = dsm(d9' reg ) = 0. If we take the divergence of 
Eq. (p]) in Landau gauge (89 = 0) we obtain, 

2n8m = Sdd = Od sm9 . (13) 

One solution of this equation is obviously, 

e ^ = -2nJ2D xy 8m y , (14) 

y 

where D xy denotes the lattice Coulomb propagator in position space. 

Instead of using the Abelian links {6 X ^} one can evaluate observables from the 
monopole and photon parts {6^ 9 } and {6^ e9 }, separately. We will refer to such expec- 
tation values as monopole and photon contributions, respectively. 

4-2. Successes 

I will present some facts in support of the superconductor confinement scenario in 
the MAP: approximate Abelian dominance of the static potential has been verified.^ 
In a recent study on a large lattice,c3 the APSU(2) string tension has been confirmed 

■fl assume the quantities to be given in lattice units a. k has dimension a -3 , F has dimension a -2 . 
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Fig. 2. Photon and monopole contributions V reg and T/ sin;? to the APSU(2) potential! V Ab 
in units a ~ 0.081 fm. 

to account for (92 ± 4)% of the SU{2) string tension. It is an open question whether 
the agreement will improve as the continuum limit is approached. 

The photon part of the potential does not contain a linear contribution. Monopole 
dominance and the factorisation V Ab {r) = V smg (r) ±_V re9 (r) approximately hold. In 
Fig. we display the result of a recent lattice study Everything is plotted in lattice 
units a ~ 0.081 fm. The monopole contribution amounts to (95 ± 1)% of the Abelian 
string tension. 

Approximate Abelian and monopole dominance has been confirmed for the light 
hadron spectrumcj of SU(3) gauge theory. In Fig. §, I display recent results on the 
nucleon and p masses that have been obtained on a lattice with extent La ~ 2.8 fm 
and spacing a w 0.175 fm. On the horizontal axis, the squared pion mass (in units 
of the string tension) is displayed, which changes as the quark mass is varied. Results 
for SU(3), APSU(3) and Abelian monopole contributions as well as for SU(3) with 
the photon part removed lie on the same curve. In the same study, it and p masses 
are found to become degenerate, as soon as the monopole contribution to the gauge 
fields is subtracted, both in SU(3) chromodynamics and the Abelian projected theory, 
i.e. Abelian monopoles are required for chiral symmetry breaking. Consistent results 
have been found by Lee et alffl and Bielefeld et al.B^ Of course, one can also argue 
that since instantons are always accompanied by Abelian monopoles,Ej removing these 
monopoles naturally results in a trivial vacuum topology. 

The one-loop /^-functions in APSU(A r c ) and SU(N C ) gauge theories agree,t3 i.e. at 
weak coupling, masses obtained in both theories should be proportional to each other. 
A recent analytical investigation has also confirmed anti-screenmgof colour fields in 
APSU(2).E3 These results are in agreement with numerical data.c3a 

For the next few observations, I refer to M. Polikarpov's talk at this conference: 
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Fig. 3. SU (3) p and nucleoli masses: Abelian, monopole and non-photon contributions.il 

monopoie^ have been found to be condensed in the confined phase of SU(2) gauge 
theory.EHc3 Furthermore, APSU(2) ppnopoles are spatially correlated with regions of 
access in the SU(2) action density^ which rules out that the monopoles are mere 
gauge fixing artefacts since their presence is reflected in gauge invariant observables. 
An effective monopole Lagrangian has been constructed which can (approximately) be 
mapped onto the Abelian Higgs model.Li 

4-3. Puzzles 

Of course, we do not expect SU (N c ) gauge theories to be identical to an Abelian 
Higgs model. So, we should become suspicious if the picture did not fail to describe the 
QCD reality at some point. Obviously, APSU(2) becomes very different from SU(2) 
gauge theory in the ultra violet. I will restrict myself to two points where something 
in the infra red might go wrong that I consider as being serious. 

The static potential between sources within the adjoint representation of SU{2) 
(adjoint potential) will saturate at a distance at which the binding energy exceeds the 
gluelump-gluelump threshold.!^ Therefore, we expect for the adjoint string tension o~ a dj 
the asymptotic value, 

r dV adj (r) 

a adj = hm — -f — = 0. (15) 
In APSU(2) we obtain for adjoint Wilson loops around a rectangle S, 

W£{S) = ~ [l + 2W Ab > 2 (S)} , W Ab ' 2 = cos \2 £ e x ,,) , (16) 



(x,/i)£dS 



where we call W Ab ' 2 the charge two Wilson loop; the adjoint source corresponds to 
a neutral and two charge two Abelian components. Obviously, as the area S goes to 
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Fig. 4. The charge two potential in APSU(2) in lattice units a ~ 0.081 fm.l 
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infinity, W A ^{S) approaches the constant value 1/3. Therefore, 
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Ab 
adj 



— lim log — ' - r 

ai^oo a W(r,t + a) 



V 



(17) 



i.e. the correct asymptotic value = is produced. 

However, lattice results indicate that the adjoint string will only break deep in 
the infra red at a distance slightly smaller than 1.5 fm,t3 a region in which we would 
expect Abelian dominance to hold if we consider the MAP theory to have relevance for 
hadronic physics. Simple models predict a linear rise of the potential in the intermediate 
region, with a slope that is proportional to the Casimir charge of the representation, i.e. 
o'adj = 8/3 cr for SU{2) gauge theory. Indeed, lattice simulations yield a linear increase, 
however with slightly smaller slope: 2a < a a dj < 8/3 a. In Fig. |4] the fundamental 



The 



potential V (r) is displayed, together with the charge two potential V 
curves are fits in accord to the parametrisation, 



AbO, i 



26 



Ab(,2) 



+ a Ab ^r. 



The upmost curve corresponds to e Ab ' 2 = 4e Ab , a Ab ' 2 = 8/3<r Ab . As in SU(2), we find 
the slope in APSU(2) to be somewhat smaller than expected. 

In conclusion, it seems that the charge two potential resembles the features of 
the adjoint potential within SU(2). However, it is not clear how we can get rid of the 
neutral contribution to the adjoint Wilson loop. For this purpose, interactions with the 
off-diagonal gluons have to be considered. 

A second puzzle that has recently been noticed^ is related to the energy dis- 
tribution within the ANO vortex between static sources. In Fig. [5], I show the 1/e 
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Fig. 5. 1/e-radius of the flux tube in SU(2) and APSU(2).il 

radius of the flux tube as a function of the source separation. The scale is set by 
a --1 / 2 « 0.45 fm. In SU(2) gauge theory only for the action density distribution a sta- 
tistically reasonable signal has been obtained while for APSU(2) we are able to display 
both, action and energy flux tube radii. While the width of the action density distribu- 
tion in SU(2) increases before it eventually saturates around a separation r « 1.2 fm 
at a value 25\/ e ~ 0.9 fm, both APSU(2) widths become consistent with a constant 
25\/ e ~ 0.36 fm as soon as r > 0.5 fm. 

Obviously, the matter fields change the structure of the flux tube and, therefore, 
affect an infra red observable. However, only the shape differs while the cross section, 
i.e. the string tension, remains (almost) invariant. The tremendous width of the SU(2) 
flux tube is somewhat counter-intuitive in view of the linear potential and the phe- 
nomenological success of string-like descriptions of hadrons; obviously, the SU{2) flux 
tube is not one-dimensional and transverse degrees of freedom should be considered. 
However, the projected theory yields what we would have expected: a thin, string-like 
object. We might speculate that the Abelian projection reveals the physically relevant 
core of the vortex that becomes obscured by the matter fields of the full theory How- 
ever, we are left with a puzzle: why does the amplitude of the underlying transverse 
string fluctuations remain that constant? One should be able to reproduce such a vortex 
in terms of an effective string theory. The simplest such theories are not renormalisable 
in four dimensions. Therefore, the width should increase logarithmically with the rate 
being controlled by an ultra violet cut-off, i.e. the scale at which longitudinal fluctua- 
tions become important and the string description breaks down. This is obviously not 
the case for our APSU(2) "string". What is the physical mechanism that prevents the 
string from widening? Is the effective string theory supersymmetricEll? What correction 
terms have to be added in order to produce such a string? 
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Fig. 7. Electric field E and magnetic super current k between two static sources. 

4-4- The Dual Superconductor in Detail 

In order to obtain an effective low-energy Lagrangian with monopoles and photons 
as fundamental degrees of freedom one can attempt to determine the free parameters 
by numerically matching the effective action to that of APSU(2)BE3 To complement 
such studies, one might probe the APSU(2) vacuum with electric (or magnetic) test 
charges to verify predictions of the effective theory and measure the values of the 
model parameters, which is the line I am going to follow here. Investigations of field 
distributions in presence of charges have been performed previously .Ea I will concentrate 
on the results from a more recent study £3 

We are probing the vacuum with static electric sources. For this purpose we 
consider three dimensional spatial cross sections (time slices) of the lattice. In Fig. || 
I have visualised where on the lattice different objects are "living" . The advantage in 
working with differential forms is that the Stokes theorem is guaranteed to be exact 
and not subject to lattice artefacts: if the differential Maxwell equations are fulfilled, 



the integrated versions automatically hold too. Of course, when finally relating lattice 
results to continuum physics, lattice artefacts re-enter the game. The (generalised) 
Maxwell equations [Eq. @] read, 

5E = p eh dE = *k- *B, (19) 
5B = Pmag , dB = *j + *E. (20) 

The charge densities p e i and p mag are the 4-components of j and k, respectively. The 
"dof -symbol denotes a temporal derivative. For a static problem, electric and magnetic 
fields decouple. In the absence of magnetic test charges, this implies p mag = 0, j = 
and, therefore, B = such that we are left with two Maxwell equations only, 5E = p e i 
and dE = *k. 

As mentioned above, one can define the monopole current via the latter equa- 
tion (Ampere's law) which is then trivially fulfilled. However, with our definition of k 
[Eq. fljjD ] neither of the equations necessarily holds. Strong fields may give rise to non- 
linear quantum corrections. Moreover, APSU(2) is not electrodynamics; in general, the 
action will be non-local and the equations of motion might differ. Our numerical study, 
however, verifies div E to disappear outside of the vicinitvof the sources while the curl 
of the electric field is identical to the magnetic current £21 I display the electric field 
for a distance r = 15a ~ 1.2 fm in Fig. [7|. The vortex is stabilised by the surrounding 
super current k. 

The starting point of our investigation is the London limit. The electric vector 
potential C on the dual lattice is defined through, 

E = *dC. (21) 

One representation of the (dual) London equations reads, 

C + A 2 k = 0, (22) 

with A being the inverse mass of the dual photon. By building the curl of this relation, 
we obtain E = — A 2 * dk which, together with the dual Ampere law results in, 

E = -X 2 5dE = A 2 (V 2 E + dp el ) . (23) 

Without an electric current j, the Maxwell equations imply, d ( B + A 2 dk) = 0, i.e. 
even in the absence of a magnetic field, the super current k remains constant and 
in general non-zero. From Eq. (|23| ) it is also obvious why A is called the penetration 
length: if we expose a superconducting probe to a homogeneous electric field, E = E z e x , 
the field strength will decay with the distance x towards the centre of the medium: 

E s (x) = E z (0)exp(-x/\)." 

We create a charge- ant icharge pair at a separation r parallel to the z-axis of our 
lattice, p e i = 5 3 (— r/2e 2 ) — 5 3 (r/2e z ). x denotes the transverse distance from the core 
of the ANO vortex. For the electric flux through any surface enclosing the charge, we 
expect $ e z = Jd 2 xE z (x, 0) = 1. For our geometry Eq. (|2~3|) reads in the centre plane 
perpendicular to the vortex, 



E z (x) = \ 2 V\E Z (x) + <M 2 (x) , E x {x) = 0. 



(24) 



A solution can be expressed in terms of the modified Bessel function K , 



E z {x) = £±K (x/\). (25) 

Within statistical accuracy, we find our data to be compatible with such a functional 
form for x > x m i n = 4.2a « 0.35 frn^J with parameter values, 

A = (1.82 ± 0.07)a « (1.3 GeV)^ 1 , $ ei = 1.44 ± 0.08. (26) 

For small x the data are overestimated by the fit since K (x) diverg — > which 
explains why the electric flux comes out to be significantly larger than one. A dual 
photon mass of 1.3 GeV is compatible with the mass of the lightest glueball in 577(2) 
gauge theory. However, the quantum numbers of a dual photon are J PC = l + ~ as 
opposed to + for this glueball. 

How can we refine our description such that the electric field not only in the 
surface region but also closer to the centre of the vortex is correctly reproduced? Obvi- 
ously, a second scale £ that is smaller than the 0.035 fm, above which the London limit 
seems to apply, has to be introduced. Such a scale appears in the Ginzburg-Landau 
(GL) equations as the coherence length of the GL wave function ^(x) that describes the 
spatial density of superconducting monopole charges, n(x) = |?/>(x)| 2 . We decompose 
if) into a phase and an amplitude /, 

^(x) = ^ 00 /(x)e ie W, /(x) x ^l. (27) 

£ characterises the decay of the monopole density towards the centre of the vortex, 
where / will vanish as superconductivity breaks down, while A controls the penetration 
of the vortex field into the surrounding vacuum. The case £ = corresponds to the 
London limit. If we increase the width of the flux tube, the dia-electric energy of the 
vortex is reduced by an amount roughly in proportion to XE 2 while the amount of 
energy we have to pay for pushing the monopoles further into the vacuum increases 
like ^E 2 . Values £ < \/2\ correspond to a negative surface energy in accord to the 
Abrikosov criterium for a classical system. This tendency to maximise the surface 
results in retardation of flux tubes: we obtain a type II superconductor while values 
K = A/£ < l/v2 correspond to a type I superconductor. From our experience we are 
prejudiced to expect a type II scenario since in electrodynamics type I flux tubes cannot 
be realised due to the absence of isolated magnetic charges. In the present situation, 
however, the presence of two isolated electric sources forces the field lines through the 
surrounding vacuum, regardless of the type of the superconductor. 

If we restrict ourselves to the perpendicular centre plane, the equation dC = *E 
implies for the azimuthal component of C (up to gauge transformations), 

Ce(x) = J- / d 2 x'E z (x') ^ pL, (28) 

ZTTX Jx'<x AlXX 
§ Along off-axis lattice directions, we obtain data for non- integer multiples of the lattice spacing. 




Fig. 8. The electric field and the amplitude of the Ginzburg-Landau wave function against 
the distance from the centre of the ANO vortex.H 



while the other components vanish. For non-constant density of magnetic charges, the 



London equation Eq. (p2|) is modified and becomes the second GL equation, 



C e {x) 



2ttx 



+ 



A 2 



-kg(x) = 0. 



(29) 



f(x)/X after having reconstructed 



We can solve this equation with respect to F(x) 
C e (x) via Eq. fl2|). 

The result is displayed in Fig. §, together with E z (x). Data obtained at x < 2.2a 
has to be treated with care since the difference between lattice and continuum versions 
of "curl" turns out to be bigger than our statistical uncertainty. For x > 4.2a the errors 
on / explode: here, no contradiction to the London limit has been found. We fit F(x) 
with the ansatz, 

F(x) = = \tanh{x/a), (30) 
A A 

which conforms to the right boundary conditions. The fit is included into the figure 
as well as the result of a fit of E z to a^ore involved four parameter ansatz that also 
respects the boundary conditions on /.e^ From the fit Eq. ( pOf ) we obtain A = 1.62(2)a. 
The fit to E z yields A = 1.84(8)a while a simultaneous fit to E z and kg yields A = 
1.99(5)a. This has to be compared with the value A = 1.82(7)a from the London limit 
fit of Eq. (pB"D. We end up with the conservative estimate, 



A 



1.84+^ a 



(0.15 ± 0.02) fm, $e, = 1.08 ± 0.02. 



(31) 



One should settle in a scaling study whether the deviation from $ e ; = 1 can be at- 
tributed to a non-trivial vacuum dielectricity constant due to anti-screening. 
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9. Effective coherence length versus distance from the centre of the ANO vertexJ 



Since the first GL equation is non-linear, we cannot consistently formulate it in 
terms of differential forms on a lattice, i.e. — unlike the Maxwell equations — we have 
to verify it in the continuum and discard data obtained at small x/a values where the 
lattice structure is still apparent. For our geometry the first GL equation reads, 



f{x) = f{x)+eh{x) 



1 2tiC 9 (x) 



l 2 



- -- x- I } fix) 
xdx \ dx 



(32) 



Our strategy is to solve this equation with respect to an effective coherence length 
£,es(%), employing the above mentioned parametrisations of E z (x) and f(x) to inter- 
polate the data. The result is visualised in Fig. ||. Results outside of the window of 
observation 2.2 < x/a < 4.2 are unreliable, for small x due to lattice artefacts and for 
large x due to lacking precision data on /. The figure contains error bands that are 
related to the uncertainties in a as well as in A. Within the window, £ c g varies by only 
10 %, i.e. the GL equation is qualitatively satisfied. Taking this variation of £ with x 



into account, we obtain the result, 



(0.251 ±0.032) fm. 



(33) 



The central value corresponds to £ = £cff(£). 

For the ratio of penetration and coherence lengths we determine, 



K 



1 

72' 



(34) 



i.e. we have evidence of a weak type I superconductor. However, we are rather close to 
the Abrikosov limit whose position might differ from 1 / a/2 due to quantum corrections. 



In order to finally settle the question of the type of the superconductor, interactions 
between two flux tubes should be investigated. This has been done in recent studies 
of the confined phase of U(l) gauge theorycj as well as in three-dimensional Z 2 gauge 
theory ,e!J the result being attraction in both cases, i.e. type I superconductivity. 

5. Summary and Open Questions 

The lattice is an ideal tool to test ideas on the confinement mechanism. Many 
infra red aspects of QCD are reproduced in the maximally Abelian projection. After 
the projection only the monopole contribution to the original gauge fields seems to 
be relevant for most low-energy properties. The dual Maxwell equations have been 
verified in APSU(2) and the fields are adequately described by the dual Ginzburg- 
Landau equations with the values A = 0.15(2) fm and £ = 0.25(3) fm for penetration 
and coherence length, respectively. These values correspond to a (dual) photon mass 
m 7 « 1.3 GeV~ 3a and a Higgs mass of mn ~ 0.8 GeV~ 2a, the ratio of which, 
k = A/£ = 0.59(13), indicates the vacuum of SU(2) gauge theory to be a (weak) type I 
superconductor. It is demanding to clarify whether flux tubes in SU (N c ) gauge theory 
as well as in the Abelian projected theory attract or repel each other. 

Electric flux tubes are found to be significantly thinner after the Abelian projec- 
tion. Contrary to SU(2) gauge theory, their width seems to saturate at a separation 
r « 0.5 fm at a value, 25i/ e ~ 0.36 fm. The Abelian projection also seems to suffer 
under problems with charges in non-fundamental representations. These observations 
require further thought. In the end we would like to clarify the role that charged gluons 
play. Ideally, one would like to completely circumvent the Abelian projection and arrive 
at a genuinely non-Abelian description of the superconductor scenario. In view of the 
existence of other reasonable proposals for the confinement mechanism more thought 
should be spend on relations between different such pictures. 

An extension of the detailed superconductor study presented from SU{2) to 
SU (3) gauge theory should be attempted. It is interesting to simulate the 4D Abelian 
Higgs model with the m 7 and m# parameters that have been predicted above and 
compare the resulting flux distributions with the result from APSU(2). This approach 
has recently been pursued by Chernodub and collaborators.^ 
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